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Despite the formal analogy with the Higgs particle, the amplitude fluctuations of the order pa¬ 
rameter in weakly-coupled superconductors do not identify a real mode with a Lorentz-invariant 
dynamics. Indeed, its resonance occurs at 2Ao, which coincides with the threshold 2 E gap for quasi¬ 
particle excitations, that spoil any relativistic dynamics. Here we investigate the fate of the Higgs 
mode in the unconventional case where 2 E gap becomes larger than 2Ao, as due to strong coupling 
or strong disorder. We show that also in this situation the amplitude fluctuations never identify 
a real mode at 2Ao, since such ’’bosonic” limit is always reached via a strong mixing with the 
phase fluctuations, which dominate the low-energy part of the spectrum. Our results have direct 
implications for the interpretation of the sub-gap optical absorption in disordered superconductors. 

PACS numbers: 74.20.-z, 74.25.Gz, 74.62.En 


In the BCS theory of superconductivity the forma¬ 
tion of Cooper pairs and their condensation into the su¬ 
perfluid state are both consequences of the spontaneous 
breaking of the 17(1) gauge symmetry at T c [ 1]. While 
the equilibrium value of the (complex) order parame¬ 
ter A 0 is responsible for a gap in the single-particle ex¬ 
citation spectrum, the fluctuations of its phase repre¬ 
sent the massless Goldstone mode that couples to the 
electromagnetic field leading to the Meissner effect. In 
addition, energetically-costly collective amplitude fluc¬ 
tuations are also possible [2-4, 4-6, 8-15]. The result¬ 
ing collective mode is sometimes named in the litera¬ 
ture ’’Schmid mode”[2, 4, 4, 6], or alternatively ’’Higgs 
mode” [8-15] to emphasise its analogy with the massive 
particle of the Standard Model[16]. 

However, such analogy holds only in part, since the 
dynamics of the superconducting (SC) Higgs mode is 
not described in general by a Lorentz-invariant (LI) rel¬ 
ativistic theory. Indeed, in conventional weakly-coupled 
superconductors the Higgs resonance occurs exactly at 
the threshold 2Ao for the proliferation of quasiparticle 
pairs that fully control its dynamics leading to a non- 
relativistic strongly overdamped mode[4, 5, 10], whose 
experimental signature emerges usually only in out-of- 
equilbrium spectroscopy[9, 15, 17, 18]. The situation 
can be different e.g. for lattice bosons at commensu¬ 
rate fillings, where fermionic quasiparticle excitations are 
absent. In this case the system has been described[6- 
8, 19, 21, 24, 25] by a SC-like LI 0(2) model: 
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where ip = (Ao + A)e l6 includes the amplitude (A) and 
phase (9) fluctuations of the SC order parameter. By 


retaining in Eq. (1) only Gaussian terms one has: 


(|A(<z)| 2 } = 


1 


, -2 


—w 2 + m 2 + c 2 q 2 


, <l%)| 2 > = 


-w 2 + c 2 q 2 ’ 


( 2 ) 

i.e. both the phase (Bogoliubov) and amplitude (Higgs) 
modes appear as elementary excitations with a LI dis¬ 
persion. The Higgs resonance, that is not overdamped[8, 
24, 25] by the decay processes into phase modes at higher 
orders, has been probed in cold atoms by properly shak¬ 
ing the optical lattice[26, 27]. Motivated by this result, it 
has been recently suggested[7, 8, 14] that the 0(2) model 
(1) is also the correct paradigm for ’’bosonic” supercon¬ 
ductors, i.e. fermionic systems where the energetic cost 
to create single-particle excitations E gap is larger than 
the SC order parameter Ao. The underlying idea is that 
the if Higgs resonance remains at 2A 0 it is untouched by 
the quasiparticle proliferation at 2 E gap , and the Higgs 
mode recovers its LI dynamics, as recently suggested for 
charge-density-wave (CDW) superconductors[10]. Apart 
from this peculiar case, in general the separation between 
Egap and A 0 occurs in homogeneous superconductors in 
the strong-SC coupling limit[28, 29]. In addition, it can 
also be induced at weak coupling by strong disorder, as 
it has been observed experimentally [30-35] in materials 
like e.g. InO^ and NbN, and understood theoretically as 
the effect of the localization of bosonic pairs with large 
7W36-38], However, in both cases it has not been 
proven yet the existence or not of a sharp Higgs mode 
at 2A 0 , which is also the prerequisite used recently[14] 
to interpret the observed sub-gap optical absorption in 
films of disordered superconductors[14, 39-41]. 

In this Letter we demonstrate the lack of a LI Higgs 
resonance at 2Ao for fermionic superconductors. At weak 
coupling, where the amplitude and phase sectors are de¬ 
coupled because of the particle-hole symmetry of the BCS 
solution, this is a well-known[4, 5, 10] effect of the quasi¬ 
particles proliferation at 2Ao, which overdamp the Higgs 
mode and spoil its relativistic dynamics. In the bosonic 
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limit, where the single-particle excitation energy E gap 
overcomes Ao, the quasiparticles are ineffective but am¬ 
plitude and phase sectors are strongly coupled, and any 
signature below E gap in the Higgs mode originates by its 
mixing to the phase mode. In the homogeneous case, 
where the regime E gap > Ao is achieved at strong cou¬ 
pling by unavoidably breaking the particle-hole symme¬ 
try, only the (Bogoliubov) sound-like phase mode sur¬ 
vives at low energy and long wavelength[2, 42-45, 47]. 
In the disordered case the mixing of the two sectors oc¬ 
curs already for weak coupling. However, the spectral 
weight found below 2 E gap in the Higgs spectrum does 
not identify a LI sharp resonance at 2Ao, questioning the 
proposed[14] relevance of the Higgs mode for the sub-gap 
optical absorption in disordered SC films. 

We start our analysis with the homogeneous case, by 
using as prototype lattice model for a superconductor the 
attractive Hubbard model on the square lattice, 


ticular at q = 0 one has x PP = —4A qF(u) and 

r \ _ 1 V"_ 4A 0 g k _ 

XpaM - N Z^E^El~{uj + i8yy (5) 

F= N ^ E^El - (w + i5Y )' (6) 


The spectrum of the collective modes is obtained as solu¬ 
tion of the equation |M(w,q)| = 0. This simplifies con¬ 
siderably in the weak-coupling limit, where x P A — 0 be¬ 
cause of the intrinsic particle-hole symmetry of the BCS 
solution[2, 5, 10, 43, 45], that enforces the integration 
over a symmetric energy interval around £ = 0 in Eq. 
(5). As a consequence the amplitude and phase sectors 
are decoupled and the matrix (S19) reduces to: 
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where U > 0 is the SC coupling strength and t the 
nearest-neighbor hopping parameter. By treating the in¬ 
teraction term in mean-field approximation the Green’s 
function in the usual basis of Nambu operators reads 
Gq 1 = iio n fjQ — £kd 3 + Ao<ti, where cq are Pauli matri¬ 
ces, A 0 is the SC gap and £k = —2t(cosk x + cos k y ) — p 
and E?k = + Aq the quasiparticle excitation in the 

normal and SC state, respectively (p being the chem¬ 
ical potential). The spectrum of the collective modes 
can be studied in the effective-action formalism]!, 3, 49], 
by decoupling the interaction term of the model (SI) 
both in the pairing and density channels by means of 
the Hubbard-Stratonovich transformation. After inte¬ 
grating out the fermions one is left with S = Smf + Sfl, 
where Smf is the mean-field saddle point action and Sfl 
the effective action for the collective degrees of freedom 
p,A,&, connected to the electron density, and the SC 
amplitude and phase, respectively[49]. By retaining only 
the Gaussian terms in the expansion of Sfl, and inte¬ 
grating out the density field p , the effective action for the 
SC degrees of freedom[49] reads Sab = Yhq V W (g)MT(g), 
where ^(g) = (A(g), 9(q)) and M at T = 0 and q —» 0 is 
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Here D s = A. £) kl/ y^r(l-Ck/A k ) is the phase stiffness, 

Xij{Q) = (7/^) E k ,i U „ Tr [Go(fc)^Go(fc - 9)^'] is the 
BCS response function with the identification by —> A 
and «j 3 -4 p, while \ PP = Xpp/) 1 + ( U/2)x PP ) is the one 
dressed at RPA level by density fluctuations [49]. In par- 


where we also used the result, valid in the BCS limit, 
F(ui —> 0) = AT f /(4A^), so that ~X PP — k = Np/( 1 — 
UNp/2), Np being the DOS at the Fermi level. From 
the M§f s = (l^l 2 )^ 1 element of the matrix (S19) one re¬ 
covers the Anderson-Bogoliubov sound-like mode, which 
is the only expected phase-density mode at T = 0[50]. 
The vanishing of the Mi^S = (|A| 2 ) -1 element identi¬ 
fies instead the pole of the amplitude fluctuations. As 
one can easily see from Eq. (6), as u —> 2A 0 is 
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so that the real part of the inverse Higgs propagator 
vanishes as i/4Ag — w 2 at the Higgs mass m = 2A 0 , 
which coincides here with the threshold for two-particle 
excitations. In addition, ImF(w) ~ (w 2 — dAg)” 1 / 2 at 
lo > 2Aq , leading to the usual overdamping of the Higgs 
spectral function in the BCS limit[4, 5, 10]. Unless one 
considers specific band structures [8], the same holds in 
general also at finite momentum, since the Higgs pole 
remains inside the quasiparticle continuum. 

The BCS result (7) shows that even when the particle- 
hole symmetry leads to a vanishing of the first-order time 
derivatives[12], i.e. the off-diagonal terms of Eq. (S19) 
and (7), the LI dynamics of the Higgs mode is prevented 
by the dynamical contribution of quasiparticles, encoded 
in the singular function F(u>). On the other hand, in 
the unconventional situation where the quasiparticle con¬ 
tinuum 2 E gap moves away from 2A 0 one expects that 
F( 2A 0 ) — const. Thus, within the diagonal structure 
(7) for the collective modes one would conclude that the 
Higgs mode should recover a LI dynamics, as recently 
suggested[7, 8, 14]. However, in the limit 2E gap 2A 0 
the dynamics of collective modes is no more described 
by Eq. (7). To show this let us first consider the ho¬ 
mogeneous strong-coupling limit U t, where the SC 
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FIG. 1: Intensity map of the spectral function of the ampli¬ 
tude and phase modes in the homogeneous case at U = 8t for 
n — 0.1 (upper panels) and n = 1 (lower panels). The solid 
red lines mark the value of 2Ao and the dashed black lines 
the value of 2 E gap . At n = 0.1 the bosonic limit E gap > Ao is 
reached, but the amplitude mode is strongly mixed with the 
phase, which dominates any sub-gap structure of the Higgs 
(with zero spectral weight at w = 0). At n = 1 the system 
is always particle-hole symmetric so the amplitude and phase 
sectors remain decoupled. However, here E gap = Ao and the 
Higgs mode has a broad spectral weight only above 2 E gap - 
The bending back of the phase mode at Q = (n,w) (where 
its weight is zero) is a signature of the degenerate CDW in¬ 
stability present in the model (SI) at n = 1 [9]. 


order parameter and chemical potential for the model 
(SI) read[9, 43, 49]: 


A o - \/n(2 — n), (9) 


At n ^ 1, as soon as the chemical potential goes below 
the band edge, E gap = a/Aq + (4 1 — \^i\) 2 becomes larger 
than Ao. However, in this limit the particle-hole symme¬ 
try is strongly violated, leading to a large XpA mixing 
between the amplitude and phase sectors in Eq. (S19). 
This can be easily understood from Eq. (5) in the limit 
t/U ~ 0, where |/r| A 0 and xpA — 4A 0 fiF(oj), with 
F(w) ~ 1/{E(4E 2 -oj 2 )) and E = V A o + /.i 2 ~ (7/2, see 
Eq. (9). In this situation the determinant of the matrix 


(S19) at q = 0 is given by 
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As one can see, the coupling XpA between the amplitude 
and phase sectors removes completely any signature at 
2A 0 and the only solution of \M\ = 0 is the (Bogoliubov- 
like) phase mode[2, 42-45, 47]. One can easily see[49] 
that this result holds also in the presence of long-range 
interactions, that modify the prefactor in the first line of 
Eq. (S22) lifting the phase mode to the plasmon, but do 
not alter the combination in square brackets, responsible 
for the disappearance of the pole at 2A 0 . By retain¬ 
ing a finite t value in the evaluation of all the response 
functions (5)-(6) one can show[49] that the spectral func¬ 
tion of the Higgs mode still preserves some small spectral 
weight above 2 E gap , that is however strongly suppressed 
with respect to the weak-coupling case. These analytical 
estimates are confirmed by the numerical computation of 
the amplitude and phase spectral functions [49] at finite 
oj and q, shown in the upper panels of Fig. 1. As one can 
see, any sub-gap feature in the Higgs spectral function 
is present only at finite q and it comes from the mix¬ 
ing to the phase, with no signature at the energy 2A 0 , 
marked by the red line. On the other hand at half filling 
(n = 1) jU = 0 so that it is always E gap = A 0 . More¬ 
over, since the particle-hole symmetry is preserved at all 
orders in U, it is XpA = 0, the amplitude and phase sec¬ 
tors remain decoupled and the Higgs mode shows only 
a weak spectral weight above 2E gap , see lower panels of 
Fig. 1. Notice that results from previous work[55] sug¬ 
gest that RPA describes reasonably well the elementary 
(Gaussian) collective excitations in the symmetry-broken 
state at strong coupling. Thus, even though one cannot 
exclude the presence of additional resonances at higher 
order, analogous e.g. to that found in the half-filled 
Bose-Hubbard model very near the superfluid-insulator 
transition[7, 21, 24], they should not be interpreted as 
signatures of an elementary amplitude mode. 

A second possible route[37, 38] to achieve the separa¬ 
tion between E gap and Ao in the model (SI) is by intro¬ 
ducing disorder as a random on-site energy V r uniformly 
distributed in the interval [—Vo, Vo]. Indeed, already at 
the level of the mean-field inhomogeneous Bogoliubov- 
de-Gennes equations [37] one sees that while the average 
order parameter (A 0 ) decreases as Vo increases, the spec¬ 
tral gap 2 Eg ap in the average density of states saturates 
to a finite value, see insets of Fig. 2, signaling the forma¬ 
tion of local boson pairs[36-38]. However, once more this 
does not imply that the Higgs mode emerges as a sharp 
resonance at 2(A 0 ), but instead it acquires sub-gap spec¬ 
tral weight due to the mixing to the phase mode, that is 
induced in the disordered case even for weak SC coupling. 
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FIG. 2: Spectral functions of the amplitude and phase modes 
at q = 0 for n = 0.875 obtained on a 20x20 lattice after 
average over 50 disorder configurations for two values of SC 
coupling (U/t) and disorder ( Vo/t ). The vertical dashed lines 
and green arrows mark the spectral gap 2 E gap and 2(Ao), re¬ 
spectively, whose disorder dependence is shown in the insets. 
For comparison we also show <r(w) (blue line) from Ref. [5]. 


the excitation of a single phase mode (phason), see Fig. 
3a, while at higher order the convoluted process between 
one phason plus one Higgs mode is also possible, see 
Fig. 3b. In the clean case the one-phason process (a) 
contributes only to the delta-like superfluid response at 
w = 0, while the process (b) contributes in general to 
<t(w) at finite frequency, with an intensity and shape that 
depends on the coupling D s of Eq. (S27) and on the spec¬ 
tral function of the Higgs mode. More specifically, in the 
clean bosonic model (1), where D s = D s = c 2 Aq and the 
Higgs mode (2) has a sharp resonance at ui = m, the pro¬ 
cess (b) leads [6-8] to a finite-frequency optical absorption 
with an edge exactly at m. By applying this result to su¬ 
perconducting films the authors of Ref. [14] interpreted 
the sub-gap optical absorption found at strong disorder 
as an indirect evidence of the existence of a sharp Higgs 
mode at an energy m < 2 E gap . However, our calcula¬ 
tions do not support this conclusion, since we explicitly 
demonstrated that for a realistic fermionic model with 
disorder the spectral function of the Higgs mode does not 
identify any sub-gap LI resonance. In addition, we also 
checked[49] that in the fermionic model (SI) is D s ~ 0 at 
weak coupling, strongly suppressing the optical process 
shown in Fig. 3b. 


In Fig. 2 we show the amplitude and phase spectral func¬ 
tions at q = 0 computed at RPA level[49] for two indica¬ 
tive values of SC coupling and disorder. As one can see, 
the Higgs mode acquires spectral weight below 2 E gap al¬ 
ready for Vo/t = 0.5 (upper panels), where E gap and Ao 
still coincide. For U/t = 5 this sub-gap feature closely 
follows the one seen in the phase sector, where it can be 
interpreted as a disorder-induced broadening of the sound 
mode, with a transfer of spectral weight from zero to fi¬ 
nite frequency. Since the mixing between the amplitude 
and phase modes vanishes as co —> 0, see Eq. (S19), this 
feature appears as a finite-energy maximum in the am¬ 
plitude sector. At larger disorder Vo/t = 3 (lower panels) 
the mixing of the amplitude with the phase is stronger 
at all energies, and also in this case the finite spectral 
weight of the Higgs mode below 2 E gap does not identify 
a well-defined resonance at the typical scale 2(A 0 ). 

Let us finally comment on the relevance of these re¬ 
sults for experiments in disordered films of superconduc¬ 
tors, where a dissipative absorption cr(oj) at w < 2 E gap 
and T T c has been recently found[14, 39-41], Since 
in this regime quasiparticle excitations are suppressed, 
the extra-absorption should arise from the SC collective 
modes, whose leading-order contribution to the current 
J can be written in general as[49, 56] 

J = —2e (d s W9 + 2£> s r ? V6») , r) = A/A 0 . (11) 

Since the optical conductivity is determined by the 
current-current correlation function, the leading-order 
optical process is proportional to (V0(a;)V0(y)), i.e. to 


(a) 

(vflwvfloo) 


(b) (**My)) 





(VtfA)Vfl(y)) 


FIG. 3: Leading-order contributions of the SC collective 
inodes to the conductivity, from the current given in Eq. 
(S27). Wavy, solid and dashed lines represent insertion of 
the electromagnetic field, the SC phase propagator and the 
SC amplitude propagators, respectively[6, 7]. 

On the other hand, in the presence of disorder the 
one-phason process of Fig. 3a contributes also to the 
finite-frequency optical response [5, 11, 58], making it 
the most plausible candidate to explain the experimen¬ 
tal findings[14, 39-41], In particular, for the fermionic 
model (SI) it has been proven[5] that the phase fluctua¬ 
tions give the largest contribution to sub-gap absorption, 
via a process equivalent to the one depicted in Fig. 3a. 
However, due to the non-trivial structure of the effective 
optical dipole of the phase mode, there is no simple cor¬ 
respondence between its spectral function and the optical 
conductivity, as one can see in Fig. 2 where we report for 
comparison a(co) from Ref. [5]. Note that the optical con¬ 
ductivity is the response to the local, i.e. screened electric 
field, and it is therefore determined by the current corre¬ 
lations functions irreducible with respect to the Coulomb 
interaction [5]. In other words, what contributes to cr(ui') 
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are the SC collective mode of the neutral case, which are 
the ones shown in Fig. 2. 

In summary, we studied the evolution of the ampli¬ 
tude (Higgs) mode in a lattice model for fermions. We 
showed that even when 2Ao goes below the threshold 
2 E gap for quasiparticle excitations the Higgs mode never 
identifies a sharp resonance at 2Ao- Indeed, the separa¬ 
tion between E gap and Ao, needed to remove the dynam¬ 
ical overscreening of quasiparticles to the Higgs mode, is 
only achieved by an explicit breaking of the particle-hole 
symmetry, which implies a strong mixing of the ampli¬ 
tude and phase fluctuations. Thus, any sub-gap feature 
in the Higgs mode is a signature of the underlying phase 
mode, and it does not resembles the sharp LI Higgs res¬ 
onance found in the bosonic model (1). These results es¬ 
tablish that the outcomes of the relativistic 0(2) model 
(1) do not describe in general the physics of fermionic su¬ 
perconductors even in the ’’bosonic” limit, as it has been 
sometimes suggested[7, 8, 14]. On the other hand, the 
relativistic nature of the Higgs mode can be recovered 
in the BCS limit when a CDW gap contributes to make 
E g ap larger than A 0 , as shown in Ref. [10]. An interesting 
open question is the possibility that a similar mechanism 
can hold also in systems like cuprate superconductors, 
where increasing experimental evidence [59-62] has been 
accumulating for the coexistence of SC and CDW order. 
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Supplemental Material 


EFFECTIVE ACTION FOR THE COLLECTIVE 
MODES: GENERAL FORMALISM 

Let us start from Eq. (3) of the manuscript, that we 
report here: 

H = ~ t Yl (4x C J> + h ' c •) - U Y1 C it C !i C 4 C iT- ( S1 ) 

(ij)o- * 

To investigate the physics of the collective modes around 
the mean field solution we follow the usual Hubbard- 
Stratonovich procedure [SI], as implemented e.g. in Refs. 
S2, S3. We then introduce in the action for the fermions 
a bosonic complex field which decouples the on¬ 

site interaction term of (SI) in the pairing channel. At 
T <T C one can choose to represent the superconducting 
(SC) fluctuations both in polar (amplitude and phase) or 
cartesian (real and imaginary parts) coordinates. In the 
former case the additional use of a Gauge transforma¬ 
tion on the fermionic operators makes the dependence 
of the effective action on the time and spatial deriva¬ 
tives of the SC phase explicit, and it is then more con¬ 
venient in the long-wavelength limit. The equivalence 
between the two approaches is guaranteed by the Ward 
identities, as discussed e.g. in Ref. [S4]. We will then 
decompose f/’A(r) = [A 0 + A i{r)\e iei ^ T \ where Aj(r) 
represent the amplitude fluctuations of ^A around the 


mean-field value A 0 and 9 its phase fluctuations, which 
appears explicitly in the action after a Gauge transfor¬ 
mation Ci —> Cje 1 ^/ 2 . The interaction term of Eq. (SI) 
can also be decoupled [S2, S3] in the particle-hole chan¬ 
nel by means of a second (real) bosonic field ip p = p 0 + p 
which couples to the electronic density Ha c i<? 

and represents the density fluctuations p of the system 
around the mean-field value p 0 . 

After the Hubbard-Stratonovich decoupling the action 
is quadratic in the fermionic fields that can then be in¬ 
tegrated out leading to the effective action for the fields 
A, 9 and p: 

^e//[A, 9, p] = Smf + Sfl[A, 6,p) , (S2) 

where Smf = ^ru~ + ^Tu ~ Trln(—G q *) is the mean field 
action, Gq 1 = icu n ao — £k<53 + Ao<ti is the BCS Green’s 
function and 

S FL . X AAS! (S3, 

z —' n 

n> 1 

is the fluctuating one, with the trace acting both in spin 
and momentum space, where E*,*,/ denotes the self-energy 
for the fluctuating fields, which reads explicitly: 


Sfefc' = ~\[^ A ( k “ k ') a i “ \f^P( k ~ k ') a 3 - \[^\ e ( k ~ “ fc, )oO-3 - (£k - CkO^o] - 

- -^J2 d ^ qi '>9(l2)^ L sin C ^sin C ^a 3 6{q 1 +q2- k + k')+ 0 (9 3 ) , (S4) 

qi,v v ~ ~ 


with k = (if2 ra ,k) and f l n = 2nTn bosonic Matsubara 
frequencies. Notice that the last line of Eq. (S4) repre¬ 
sents the transcription on the lattice of the usual (Vd ) 2 
term for a continuum model. In addition, in contrast 
to the continuum model, the lattice self-energy (S4) de¬ 
pends in principle[S2, S3] on all higher-order powers of 
the 9 held, which are however irrelevant for the derivation 
of the Gaussian action. 

To derive the Gaussian action for SC fluctuations we 
should retain the terms up to n = 2 in Eq. (S3). The 
terms coming from an insertion of the a 0 term of Eq. 
(S4) describe the effects of a paramagnetic current, so 
that they lead for example to the depletion of the super¬ 


fluid stiffness at finite temperature due to quasiparticle 
excitations. On the other hand, mixed terms containing 
a do times a <ti, <73 matrix give higher-order contribution 
in q. Thus, since we are interested in the T = 0 and long- 
wavelength limit, in the clean case we can safely neglect 
these terms. With lengthy but straightforward calcula¬ 
tions one can then show that at gaussian level Sfl can 
be written as: 


S FL [A,9,p] = ±'52&(q)M FL (q)*(q) - (S5) 

q 
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with 4 t T (q) = (A (q) 6(q) p(q )) and: 


with M given by Eq. (4) of the manuscript, i.e.: 


/2/U + \AA(q) l ~fXpA{q) X P A (q) 

Mfl=\ -fxpA(-g) tXw.(«) + i 1 ®(q) ^¥Xpp(9) 

V x P A(-q) l fx P p{q) 2/u + x P p(q) 

(S 6 ) 

Here Xij (<?) = w, Sfc Tr [ G o(k + < 7 )cr, ; G 0 (fc)<T.,] are the re¬ 
sponse functions computed at BCS level and evaluated 
in the zero temperature limit. Since in Eq. (S4) the in¬ 
sertion of a < 7 i or <73 Pauli matrix corresponds to a term 
proportional to the amplitude or to the density/phase 
fluctuations, respectively, we made the correspondence 
1 — > A, 3 — > p in the notation for the BCS susceptibil¬ 
ities that appear as coefficients of the action (S5)-(S6). 
They are explicitly given at T = 0 by: 


Xaa (q) 


1 y^ E k + E k . E k E k + ~ A o 

N ' k (w + iS) 2 - (E+ + E^)' 2 ’ 

(S7a) 


X P A{q) 


A 0 y^ E w + E k ._ £k + £k _ 

N 'k E i E k ‘ (w + i<5) 2 -(£;+ +s-) 2 ’ 

(S7b) 


/ , J_ y^ -^k + ^k -^k-^k ~ ^k^k + A o 

" Ns ^ Ke - ■ {uJ + i6)2 _ {Ei + E -^ 

( S7c) 

with = i?k± q /2 an d £k — Ck±q/2- 

The quantity D s (see also Eq. (S29) below) is the 
phase stiffness, that appears as coefficient of the term 
to(q) = 4 sin 2 (g l/ /2) ~ q 2 , so that in the hydro- 

dynamic limit the inverse of the bare phase-fluctuation 
propagator is: 4M 22 (g r ) ~ —k< 7 q + d? s q 2 , which de¬ 
fines the BCS sound velocity: = \JD s /k , where 

k = —Xpp(O) = i^^k-^k the bare compressibil¬ 
ity. By means of the self-consistent equation for Ao, i.e. 
2/U = ]T) k l/i?k one can also rewrite the M p l term as: 


2 /U + ,\aa(w, 0) 


1 w 2 — 4Ag 

N~s^E k [(w + isy - 4El] = 

(4Aq — w 2 )F(u) (S8) 


that corresponds to the expression used in Eq. (4) of the 
manuscript. To describe the fluctuations only in the SC 
sector we can integrate out the density fluctuations p in 
Eq. (S 6 ), so that the Gaussian action for SC fluctuations 
reads: 

Sfl[ A, 6\ = 1 £ (A(-g) 0{-q)) M(q) , 

(S9) 


M(q) 


(2/U + Xaa(q) ^ tfXpA(q) \ 
V - l fx P A(-q) ^-x PP (q ) + T w (q)/ 


t (S10) 

The integration of the density field is equivalent as usual 
to the RPA dressing of the BCS susceptibilities. In par¬ 
ticular, one has that x.ab = Xab — us i n S the 

equivalence 1 /(2/C/ + xpp) — Xpp/Xpp one is then left 
with Eq. (4) of the manuscript. 


EXPANSION AROUND 2 E gap AT STRONG 
COUPLING 


Let us investigate in detail the amplitude fluctuations 
at q = 0. From Eq. (S10) we see that at q = 0 the 
RPA resummation of the bubbles in the density channel 
factorizes out, so that the spectral function of the Higgs 
Pa(w) = dlm{l/XAA(w + id, q = 0 )} is determined by 
the frequency behavior of the function: 

Aaa = -^+Xaa-^. (Sll) 

U Xpp 

Let us then investigate its behaviour from weak to strong 
coupling. As we discuss in the main text, at strong cou¬ 
pling one can reach the bosonic limit E gap ~ U/2 A 0 
only away from half filling, by moving the chemical po¬ 
tential n ~ —17(1 — n)/2 below the band edge. In the 
following we will consider for instance the case where 
p < 0 goes below the band edge —At at strong cou¬ 
pling. Thus, by introducing E m ax = ^A^ + (At + \p\) 2 
and E min = y/A% + ( 4t ~ M ) 2 one has that at strong 
coupling Eg a p = E m i n . If we also approximate the mo¬ 
mentum integration in the Eqs. (S7) with an energy in¬ 
tegration over a constant density of stated Np we can 
put: 


^+Xaa = (4A 2 - w 2 )/ 0 (w), (S12a) 

XppM = 4Aq/ 0 (w), (S12b) 

Xpa(w) = 4A 0 /i(w), (S12c) 

where the real and imaginary parts of the Io,i functions 
are given explicitly by: 
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ReJo(w) = — V 




Np 


J-H-At V^O 

AT f 0(2A o — to) 


V / A2TF( W 2 - 4Aq - 4 £ 2 


+ 


oi\J 4Aq — w 2 

JV>0(w - 2 Aq) ^ 


arctan 


( (4* + \fi\)u 


- arctan 


l {4t- \n\)u> 
V V 4Aq - u 2 E 


V ^4A 2 - oj 2 E, 

(4 1 + \]j\)u_ + E max sJu> 2 - 4A(j (4t - |^|)o; + E min sJ w 1 - 4A^ 


-/ — in -. — ' - . — 

- 4A(5 (4( + - A',,.,, , y o.'- - I A/, (4( - A - v .. - - IA,: 

f W<4(, 

\ e (“- 2£ ~)i^fc3si' l"l> 4 ' 


r» —//+4£ 


= -p 


Re/i(w) 

Ini/, (A) — @(cc 2 P/ m?n ) 




-M-4t v / Ag+F(^ 2 -4A2-4e 2 ) 


Np 

4 w 


In 


CU H~ 2 Ey, 


ui — 2 EL 


w — 2£L 


w ■ 


■ 2EL 


4w 


> 4< 


(S13a) 

(S13b) 

(S13c) 

(S13d) 

(S14a) 

(S14b) 


so that e.g. Io is an approximated expression for the 
function F(co) introduced in Eq. (S 8 ) above. 

In the weak-coupling regime the chemical potential 
lies in general well inside the band \p\ <C 4t, so that 
E max ~ E m i n and the integral (S14a) is approximately 
zero. This leads to the well-known decoupling of the 
amplitude and phase sectors due to the (approximate) 
particle-hole symmetry of the BCS solution, already 
mentioned in the main text. In addition, as far as 
|/i| < 4t the two arctan in Eq. (S13b) have the same 
sign as ui —> 2 A(/, so that they give a constant contri¬ 
bution and Reio(w) diverges as 1/-\/4Aq — w 2 , leading 
to ReXAA(w) ~ (w 2 — 4Ag)/ 0 (w) oc ^/4Aq — w 2 at the 
Higgs pole. At the same time Re / 0 is finite asw-> 2Aq 
since the logarithmic term in Eq. (S13c) vanishes ex¬ 
actly as i/w 2 — 4A 2 , so that RcAaa(w) ~ (w 2 — 4Aq) 
as oj > 2Ao- On the other hand when oj > 2Ao the 
imaginary part Iml 0 (u;) diverg es, see Eq. (S13d), so that 
IiuA'aa grows as \J oj 2 — 4Aq as oj > 2 A 0 , leading to 
the typical over-damped resonance at 2A 0 in the spec¬ 
tral function p/\, as shown in Fig. Sla. It is worth noting 
that even retaining the small but finite value of Re/i(w) 
does not change these results. Indeed, as far as p lies in¬ 
side the band the function I\{oj) (and then x p a) remains 
finite around oj = 2A 0 , see Eq. (S14a), while xpp is P r0_ 
portional to Io, leading to the same results discussed so 
far in the case of perfect particle-hole symmetry. 

At strong coupling (and away from half filling) as soon 
as |p| > At the optical gap moves from 2A 0 to the higher 
value 2E gap = 2 E m i n . In this case the two arctan in Eq. 
(S13b) have opposite sign as oj —> 2A)/, so that they give 
overall a contribution proportional to a/4Aq — oj 2 that 
removes the divergence of Re/o(w = 2A 0 ). At the same 
time Im / 0 is only finite above 2 E gap , see Eq. (S13d), so 


that the amplitude fluctuations described at BCS level by 
the single (S12a) term would recover a perfect Lorentz- 
invariant (LI) dynamics, with a sharp spectral function at 
2Ao- However, as soon as p moves outside the band edge 
the particle-hole symmetry is strongly violated and the 
coupling to the phase, dictated by the XpA oc / 1 function, 
becomes large, moving the spectral weight of the Higgs 
fluctuations away from 2A 0 towards the new optical gap 
2E gap . In particular one can see from Eqs. (S13c) and 
(S14a) that both Re/o and Re/i diverge logarithmically 
as oj —y 2 E m i n = 2E gap , so that! 


Re/oM 


N F 

8 E gap (\p\ - 41 ) 


In 1 - 


— ) 

4 E' 2 / 

gap / 


+ Kq ; 

(S15) 


ReE^w) 


Np 

QEg a p 


In 1 - 


— ) 

4 E 2 / 

gap / 


+ A'i 


(S16) 


where Kq = 
and Ki = 


8 Eg, 
1 


H ln 

I E, 


4(//—4 1) 2 (fi+4t)E gap — (n~4t)E„ 
— E a 


8 Eq 


In ^4 _ According to Eqs. 


(S15)-(S16) now all the bubbles entering the definition 
of the inverse Higgs propagator (Sll) are singular at 
2 Eg a p, so that at oj < 2 E gap : 


Xaa — 


-4A 2 -4 ± 


~ 8(p — 4t) 2 K 0 + 8(p — 4t)K\ 


(S17) 


where we used the fact that w 2 — 4Ag — 4(/i // 0 ) 2 ~ 
w 2 — 4 Eg ap (plus 1/| ln(. ,.)| terms) so that the diver¬ 
gence of /o(w) as u —> 2E gap is compensated by the 
quantity in square brackets in Eq. (S17), and only a finite 
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<W(2A 0 ) 
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FIG. SI: (a) Frequency dependence of the spectral function 
for the amplitude fluctuations at q = 0 at U = 8t for increas¬ 
ing value of the chemical potential. The numerical results 
are obtained with a broadening <5 = O.lAo. The scaling fac¬ 
tor No represents the equivalent density of states needed to 
match the numerical result with the analytical expressions 
discussed in the text, and it accounts for the differences due 
to the change of doping. The inset shows the evolution of the 
SC order parameter Ao. For comparison we show in panel 
(2) the spectral function of the BCS term (S12a) alone, that 
neglects the mixing to the phase sector. As one can see, as 
soon as the chemical potential moves outside the band edge 
(|/r| > 4t) the BCS spectral function in panel (b) displays 
a strong sharpening, since the quasiparticle contribution en¬ 
coded in the function Jo(w) is no more diverging at w = 2Ao, 
and the (bare) Higgs mode recovers a LI dynamics. How¬ 
ever, as shown in panel (a), the coupling to the phase mode 
removes completely the signature at 2Ao in the full Higgs 
spectral function and only a weak signature at 2 E gap > 2Ao 
is left. The inset of panel (b) is a zoom around 2 Aq. 


value remains. This result coincides with the simplified 
expression computed at t = 0 and quoted in Eq. (10) 
of the main manuscript. The finite value of ReA'AA as 
oj —> 2 E gap , and the finite values of the imaginary parts 
of the / 0 , 1 1 functions at 2 E gap , imply an even weaker res¬ 
onance of the Higgs spectral function at the quasiparticle 
threshold 2 E gap when compared to the BCS case. This 
is clearly shown in Fig. SI, where we report the spectral 
function of the Higgs pa(w) computed at q = 0 with a 
numerical computation of the full expression (Sll) on the 
lattice model. As one can see, the Higgs spectral func¬ 
tion is always overdamped and, as |/i| exceeds the value 


41, its maximum moves away from 2Ao towards slightly 
higher frequencies, with a further weakening and broad¬ 
ening of the optical-gap signature. For comparison we 
also show in Fig. Sl(b) the BCS spectral function, i.e. 
the bare term (S12a) in the inverse Higgs propagator, 
pg cs (w) = ilm{l/(2/[/ + Xaa(w + iS, q = 0))}. Here 
as soon as p goes below the band edge the quasiparticle 
continuum moves away from the Higgs pole leading to 
a sharp LI resonance. However, the unavoidable mixing 
to the phase removes completely this signature, leading 
back to the broad spectral function shown in panel (a). 


ROLE OF LONG-RANGE INTERACTIONS 


The basic mechanism discussed in the manuscript (see 
also previous Section) leading to to the disappearance of 
the Higgs pole at 2A 0 in the strong-coupling case holds 
also in the presence of long-range interactions, which 
modify however the nature of the phase mode. The ef¬ 
fect of Coulomb interactions can be taken into account 
by adding to the Hamiltonian (SI) an interacting term 

H° = 2 ^( t l) c k+q,o- C k'- q ,CT' C k' > CT' c k ) o- (S 18 ) 

k,k / ,q 


where E(q) is the Fourier transform of the Coulomb po¬ 
tential in the D dimensional lattice. At small q it re¬ 
duces to the expression in the continuum limit, so that 
R(q) —> Ae 2 /|q| £)_1 where A = Air/es for D = 3 while 
A = 2tt/€b for D = 2, eb being the background dielectric 
constant. This term can be decoupled by means of the 
same HS field ip p = p a + p introduced in Eq. (S2) to de¬ 
couple the Hubbard term of Eq. (SI) in the particle-hole 
channel. As a consequence, one can follow the same steps 
which lead to Eq. (S10) with the only difference that af¬ 
ter integration of the density field the BCS bubbles are 
RPA dressed with the overall potential U — 2C(q). The 
matrix M matrix describing the coupled amplitude and 
phase fluctuations is then written as Eq. (S10), corre¬ 
sponding to Eq. (4) of the main text: 


M = 


^(4Aq - w 2 )em ^ 


iuj XpAXpp 
2 


iuj XpAXpp 
2 Xpp 


Xpp 


uj -C/LR 
-rXpp 




where the dressed charge susceptibility x pp reads: 


y lr = 

Xpp — 


Xpp 


1 + (t -v(q)) xpp 


(519) 

(520) 


If we now consider the long-wavelength limit we see that 
X RR —t —1/V r (q) 0, so that the matrix (S19) reduces 
to: 


M = 


(4Aq — w 2 )F(ui) 

0 


Xpp 


4V(q) 



(S21) 
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In the weak-coupling BCS limit where XpA — 0 one then 
recovers the equivalent of Eq. (7) of the main text, with 
the difference that now the phase mode is lifted to the 
plasmon. In particular, by considering the q —► 0 limit 
of the Coulomb potential one then finds that the plas¬ 
mon dispersion is u>p(q) — \J ( 2'Ke 2 D s /e B )q in D = 2 
and cop(q) — yj ^Tre 2 /e B in D = 3. On the other hand in 
the strong-coupling regime, where also \pA becomes sin¬ 
gular at 2 E gap , the pole of the amplitude mode is shifted 
away from 2Ao by the charge-amplitude coupling, and 
one recovers for \M\ the equivalent of Eq. (10) of the 
manuscript: 


\M\ = 


D 


_4V(q) 


s q 2 


(4A g- w 2 )F( W )-^ 
X PP 


4 [w 2 -Wp(g)] F(w)[4(Aq + p. 2 ) — w 2 ] ~ 

[u 2 -u, 2 P (q)\ . (S22) 


\e 2 /e B 
2 q 2 

UXe 2 /eb 


As a consequence, the Coulomb interaction only modifies 
the nature of the phase mode, but it does not change the 
mechanism responsible for the lack of a Higgs signature 
at 2A 0 in the case where E gap > A 0 . 

These analytical estimates are confirmed by numerical 
computation of the amplitude and phase spectral func¬ 
tions shown in Fig.s S2(a-b). Even though our model 
is in D = 2 we nonetheless simulate also the case of a 
3 D plasmon for the sake of completeness. We adopt the 
following form of the long-range Coulomb potential on a 
D-dimensional lattice 


F C (q) 


_ Vo _ 

[2 - cos(q x ) - cos(g y )] (£l ~ 1)/2 


which at small momenta reproduces the 1 /q 2 (1 /q) be¬ 
havior in D = 3 {D = 2) dimensions. In the 2D case, 
see Fig. S2a, the ujp ~ yjq dispersion of the plasmon 
found in the phase sector reflects also in the amplitude 
one at finite momentum. In the 3 D case, see Fig. S2b, 
the plasma mode can be pushed even above E gap for suf¬ 
ficiently strong Coulomb potential, so that no spectral 
weight appears below 2 E gap . However, in both cases one 
clearly see that no signature is left at 2Ao, in full agree¬ 
ment with the conclusions of the manuscript, that holds 
irrespectively on the character (short or long range) of 
the density interactions. 

For what concerns the inhomogeneous case one expects 
similar results: the presence of the Coulomb interactions 
does not eliminate the mixing between the amplitude and 
phase sector, but it affects in general the nature of the 
phase mode. On the other hand, since we are inter¬ 
ested in making a comparison between the Higgs spec¬ 
tral function and the optical response, see Fig. 2 of the 
manuscript, we can in first approximation omit the RPA 
resummation of the Coulomb potential. Indeed, as we 
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b) U/t=8, n=0.1,3D Coulomb potential 
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FIG. S2: Intensity map of the spectral function of the am¬ 
plitude and phase modes in the homogeneous case at U = 8t 
for n = 0.1 in the presence of a 2D (a) or 3 D (b) Coulomb 
potential (Vo /t = 100). The solid red lines mark the value of 
2 Ao and the dashed blue lines the value of 2 E gap - 


already discussed in Ref. [S5], the optical conductivity 
is the response to the local, i.e. screened , electric field 
and therefore it is determined by the irreducible current 
correlations with respect to the Coulomb potential. In 
other words, in disordered systems the Coulomb interac¬ 
tions affect the phase mode, but they do not affect its 
projection on the sub-gap optical absorption, whose un¬ 
derstanding is one of the main experimental motivations 
of our work. 


SPECTRAL FUNCTION IN THE DISORDERED 
CASE 

In the Sec. I above we outlined the derivation of the 
effective action for the SC degrees of freedom by using 
a polar-coordinates representation for the HS field ipA 
further supplemented by a Gauge transformation. In¬ 
deed, this approach makes more transparent the struc¬ 
ture of the collective modes in the long-wavelength limit, 
discussed in details in the Sec. II. On the other hand, 
for the computation of the finite-q amplitude and phase 
spectral function, and in particular for the evaluation of 








































them in the presence of disorder, it is more convenient to 
introduce the pairing operators in real space: 
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Srji = - (cqCj t ), (S23) 


and to use a cartesian representation of local amplitude 
(A) and phase (0) fluctuation operators 

SA t = ( Sr/i + Sril)/V2 
5<di = i(Sr]i - 5rjj)/V2 . 

Together with the local charge fluctuation 

Spi = ^ ' ^ C Lr C io- ~ ( C io- C icr)^ j 
a 

one can set up a matrix of correlation functions 

X^M = * J dte iut (Td n (t)R m (0)) (S24) 



FIG. S3: Dependence on the SC coupling strength U/t of the 
ratio Ds/D s appearing in the definition (S27) of the current 
for the homogeneous model (SI) at various densities. 


which can be computed on the RPA level as described in 
Appendix A of Ref. [S5]. 

For the homogeneous system one can compute explic¬ 
itly the correlation functions in momentum space. With 
respect to the bare susceptibilities defined in Eqs. (S7) 
the only difference is in the one coupling the amplitude 
and phase sector, that reads: 


Xo 


Ae 




2 N, 

1 


ji t_ 

lK 


£k , £ 


E~ 


1 


uj + iS - E k - E k 


+ iS + E k + E k 


(S25) 


As one can see, as q —>• one recovers Xo° = (w/Ao)x p a, 
see Eq. (S7b). On the other hand at large q to recover 
Eq. (S25) one should add in the expansion (S3) also the 
terms coming from the E ~ V0 <jo in Eq. (S4), negligi¬ 
ble in the long-wavelength limit. Apart from this, the 
RPA resummation (as discussed in Appendix A of Ref. 
S5) leads to the same scheme outlined in Sec. I. Indeed, 
at RPA level the present formulation in cartesian coordi¬ 
nates is completely equivalent to the polar representation 
of SC fluctuations used in the previous sections, as one 
can show explicitly by means of the generalised Ward 
identities (cf. e.g. Ref. [S4]). The imaginary parts of 
the full amplitude Xq' A ( u ) and phase Xq' 0 ( w ) correla¬ 
tion functions are the quantities shown in Fig. 1 of the 
main paper. 

To investigate the effect of disorder we add to Eq. SI 
a random onsite potential, 


H dia = E^ c Ua (S26) 

i,cr 


with — Vo < Vi < Vo being a random variable which is 
taken from a flat and normalized distribution. The cor¬ 
relation functions are then calculated on finite lattices 


(typically 20 x 20) and averaged over disorder configura¬ 
tions. The same approach can be used to calculate the 
optical conductivity in the presence of disorder (cf. Ap¬ 
pendix A of Ref. [S5]) and the corresponding results are 
shown in Fig. 2 of the main manuscript. 


SUPERCONDUCTING COLLECTIVE-MODES 
CONTRIBUTION TO THE CURRENT 


Let us outline briefly the derivation of Eq. (11) of the 
main manuscript. In the 0(2) model the current is easily 
obtained by means of the minimal-coupling substitution 
—*V —> —tV + 2eA as 


3 S = -2ec 2 (V’V-!/>* - V*VV>) = 

~ —2ec 2 Ap(V0 + 2rnS70) = 

= -2e(D s V6> + 2D S 7/V(9) (S27) 

where D s = D s = c 2 Aq and we put ip = A 0 (l + r))e l9 . 
Here we retained only leading-order terms in the am¬ 
plitude (if) and phase fluctuations, in analogy with the 
discussion of Refs. [S6-S8]. The structure of the current 
in the lattice model (SI) differs in general from the ex¬ 
pression (S27). Indeed, apart from the presence of the 
quasiparticle contribution to the current, absent in the 
bosonic 0(2) model, the collective-modes contribution to 
the current does not display in general the local structure 
encoded in the form (S27). Nonetheless, since in general 
J s = — dSFL/dA, a structure similar to Eq. (S27) can 
be derived by considering in the action (S3) all the terms 
of the form: 


Sfl ^ 8 


dx 


D s + 2D s r,(x) |V0(®)| 


(S28) 
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where D s is the superfluid stiffness, defined as usual as: 


D s = 



(S29) 


the Dg A vertex of Eq. (S30), after integrating out the 
density field at RPA level. This procedure leads indeed 
to the form (S28) of the effective action, with the effective 
coefficient 


Let us then consider again the expansion (S3) and 

(3) 

let us derive all the terms of third order in S y F p = 
Tr [GoS^GoEW] containing (V0) 2 times an other fluc¬ 
tuating field, that can be either the amplitude or the 
phase. By considering terms at second order n = 2 
in the action (S3) arising from the convolution between 
y;( 2 ) ~ 0 2 ct 3 and E^ 1 ) ~ Ady + pa 3 we have: 


S fl = \j dxdy\\79(x)\ 2 (r](y) + p(y))R(x — y) = 


dx\X79{x)\ 2 [2D eAV (x) + 2 Dg pP {x)} (S30) 


where 


Dg a 


Dg p 


^0 V' <5 2 £k £k _ Aq Y" e k?k /qoi\ 
W s ^ dki El 2 N s ^ ^ 

^0 1 _ Ap £ k 

W s ^dklEl 2N s ^El 


As one can see, in the continuum limit where ~ 
k 2 /2r?z — p one has that D s ~ n/m , Dg A — 
-(A 0 /2m)x P A(w = 0), and Dg p ~ -(A 0 /2m)x P p(w = 
0), where x P A and xpp are given by Eqs. (S12c) and 
(S12b), respectively. As a consequence, the approximate 
particle-hole symmetry of the BCS solution also guaran¬ 
tees that at weak coupling the phase gradient decouples 
from the amplitude fluctuations, while Dg p is finite but 
usually smaller than D s , since Dg p ~ (A 0 /in)Np ~ A 0 
while D s is of order of the Fermi energy. In the oppo¬ 
site high-density limit, i.e. near half-filling where p ~ 0, 
Dg p ~ 0 while Dg A is finite, but always much smaller 
than D s since Dg A ~ Ao(Ao/t) while D s ~ t, with 
Ao /t <C 1 at weak coupling. Thus, in the weak-coupling 
regime all the processes involving the excitation of one 
phason plus one Higgs or density mode are suppressed in 
the fermionic model (SI). 

In the strong-coupling regime, where x P A becomes as 
large as \pp away from half-filling, we expect that the co¬ 
efficients (S31)-(S32) become as large as D s , and a form 
similar to Eq. (S27) can be recovered. To account also for 
the contribution of the density field, absent in the 0(2) 
bosonic model, we will make the very rough approxima¬ 
tion to consider its effect only in the renormalisation of 


Ds = Dg A — 


A?pXpa(0) 


(S33) 


2/u+xpp^y 

In the strong-coupling regime the quantity (S33) can be 
easily estimated at leading order in the small coupling 
a = (2 t/U), see e.g. Ref. [S9]. In particular one has 
that: 


XpA - 


(S34) 

Xpp — 


(S35) 

D s ~ 

Ua 2 (l — S 2 ) 

(S36) 

Dg A — 

~ja 2 (l-S 2 )(l-3S 2 ) 

(S37) 

Dg p — 

^a 2 S(l - S 2 ) 3 D 

(S38) 


where S = 1 — n and Eq. (S34) only holds for d / 0. 
As a consequence one immediately sees that D s /D s —> 1 
for any finite doping away from half-filling, where the 
bosonic limit is recovered. On the other hand at half¬ 
filling xpA = 0 at all values of U so that D s = D SA 
always and D s /D —► —1/2 at strong coupling, as shown 
in Fig. S3. 
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